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ABSTRACT 

Mercury is the target of two space missions: MESSENGER (NASA) which orbit 
insertion is planned for March 2011, and ESA/JAXA BepiColombo, that should be 
launched in 2014. Their instruments will observe the surface of the planet with a high 
accuracy (about 1 arcsec for BepiColombo), what motivates studying its rotation. 
Mercury is assumed to be composed of a rigid mantle and an at least partially molten 
core. We here study the influence of the core-mantle interactions on the rotation 
perturbed by the solar gravitational interaction, by modeling the core as an ellipsoidal 
cavity filled with inviscid fluid of constant uniform density and vorticity. We use both 
analytical (Lie transforms) and numerical tools to study this rotation, with different 
shapes of the core. We express in particular the proper frequencies of the system, 
because they characterize the response of Mercury to the different solicitations, due 
to the orbital motion of Mercury around the Sun. We show that, contrary to its size, 
the shape of the core cannot be determined from observations of either longitudinal 
or polar motions. However, we highlight the strong influence of a resonance between 
the proper frequency of the core and the spin of Mercury that raises the velocity field 
inside the core. We show that the key parameter is the polar flattening of the core. 
This effect cannot be directly derived from observations of the surface of Mercury, but 
we cannot exclude the possibility of an indirect detection by measuring the magnetic 
field. 

Key words: planets and satellites: individual: Mercury - planets and satellites: 
interior 



1 INTRODUCTION 



Mercury is the target of two current space missions (see e.g. iMcNutt et alj (|2004f) ). The first one, MESSENGER (NASA), 
already performed three flybys on January 14, October 6, 2008, and September 29, 2009, before orbit insertion in March 
2011. The second one, BepiColombo (ESA/JAXA), is planned to be launched in 2014 and to reach Mercury in 2020. The 
preparation of these two missions motivated an in-depth study of the rotation of Mercury. 

The rotation of Mercury is a unique case in the Solar System because of its 3 : 2 spin-orbit resonance, Mercury performing 
exactly 3 rotati ons during 2 revolutions about the Sun jPettengill and Dvcdll965h . It corr esponds to an equil ibrium state 
1 Colombd 119651 ) known as Cassini State 1. Recently, radar Earth-based measurements by iMargot et al.l (120071 ) detected a 
88-day longitudinal libration of Mercury with an amplitude <j> of 35.8 ± 2 arcsec. This amplitude being nearly twice too high 
to be consistent with a rigid Mercury, it is the signature of an at least partially molten core. If we consider Mercury as a 
2-layered body with a rigid mantle and a spherical liquid core that does not follow the short-period (~ 88 days) excitations 
and does not interact with the mantle, we can derive from this amplitude the inertia of the mantle plus crust. In part icular, 
nami ng C' m the inertial polar momentum of the mantle and A < B < C the inertial momenta of Mercury, we have ijPeale 
1972h : 
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Figure 1. In the left panel we have 3 reference frames: one linked to the ecliptic plane (el, ej, e$), another linked to the angular momentum 
N (ni , , njj ) , and the last one linked to Mercury's axes of inertia {fi,f'2,f:j)- In the right panel we have a similar configuration but 
instead of the angular momentum of Mercury, we have a reference frame linked to the angular momentum of a pseudo-core (defined 
later). We have the Euler angles (h,K,g) positioning the vector 712 on the plane perpendicular to the angular momentum of Mercury 
and the Euler angles (h c , K c ,g c ) positioning the vector n!j on the plane perpendicular to the angular momentum of the pseudo-core. The 
angles (I, J) and (l c , Jc) position the axis of least inertia. Note that J c is defined on the other side than J. 
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where C22 is a second-degree coefficient of the gravitational potential of Mercury, M its mass, R its radi us, and e its orbita l 
eccentricity. It leads to (B - A)/C m » (2.033 ± 0.11 4) x 10" 4 with e » 0.206. If we take C/MR 2 = 0.34 jMilani et alj|200lh 
and C22 = (1.0 ±0.5) x 10 ~ 5 jAnderson et alj[l987h we get C m /C = 0.5 79^ ^. 



Recent studies in one l|Peale et al]|2009l ) and two jPufev et al-lboosh degrees of freedom have theoretically estimated the 



longitudinal librations of Mercury. They highlighted in particular the possibility of a resonance with the jovian perturbation, 
whose period is 11.86 years, that could potentially raise the amplitude of a long-term (~ 12 years) libration. Other periodic 
terms of a few arcsec have been es t imated. This model also predic t s that the latitudinal motion of Mercury should adiabatically 
follow the Cassini State 1 ( Pealei 1 20061 ). D'Hoedt and Lemaitre! ( 20081 )). with short-period librations of about 10 milli-arcsec 
1 Dufev et al.ll2009). In all these stud ies, the core-mantle interactions are neglected. 

Recently. iRambaux et all J2007I) explor ed the dynamics of the rotation of Mercury, including core-mantle interactions in 
t he SONyR model llRambaux and Bolsll2004l) . We here propose an alternative study, starting from the Hamiltonian formulation 
of lTouma and Wisdom! ((2001) and highlighting the dynamical implications of core-mantle interactions, by considering Mercury 
as composed of a rigid mantle and a triaxial ellipsoidal cavity filled with inviscid fluid of constant uniform density and vorticity. 



2 THE INTERIOR MODEL 



The differe ntial equation s ruli ng the motion o f a 2-layered bo dy with a rigid mantle and a liquid non-spherical core have been 
derived by llloughl J 18951) andlPoincarel jl910l ). More recently, iTouma and Wisdom! (|200ll ) gave a Hamiltonian formulation of 
this problem, that lHenrard (2008) applied to the rotational dynamics of Io, assuming that the core and the mantle were aligned 
and proportional. Here, we generalize the model of Henrard, allowing the core to be non-proportional and non-spherical. 



2.1 Physical model 

Four references frames are being considered (see Fig[T]&[2]). The first one, (e~i, e~2, e|) is assumed to be inertial for the rotational 
dynamics, it is in fact centered on Mercury and in translation with the inertial reference frame in which the orbital ephemerides 
of Mercury are given. This reference frame is related to the ecliptic at J2000. The second one, (nf,n2,n§) is linked to the 
angular momentum of a pseudo-core that we define later, while the third one, i.e. (ril, n®, ^3), is linked to the total angular 
momentum of Mercury. Finally, the last one, written as (/1, /b, fa), is rigidly linked to the principal axes of inertia of Mercury. 
In this last reference frame, the matrix of inertia of Mercury reads: 




(2) 



with < A ^ B ^ C, while this of the core is: 
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in the same reference frame. So, the orientations of the mantle and the cavity are the same, a misalignment of their principal 
axes would require to consider the mantle as elastic, this is beyond the scope of the paper. As for the whole Mercury, we have 
< A c ^ B c ^ C c . In this way, the principal moments of inertia of the mantle are respectively A m — A — A c , B m = B — B c 
and Cm — C — Cc The principal elliptical radii of the cavity are written respectively a, b, c, yielding 



Ac 
B c 

Cc 



{x\ + x\\)pdxi dx2 dx3 
{x\ + x\)pdxi dx2 dx3 
Jj (xi + x\)pdx\ dx2 dx 3 



(a 2 + c 2 ), 
(a 2 + 6 2 ), 



where p is the density of mass of the fluid core, the integration being performed over the volume of the core. 



2.2 The kinetic energy of the system 



A Hamiltonian formulation of such a problem is usually composed of a kinetic energy and a disturbing potential, here the 
solar perturbation. Therefore, we consider every internal proc ess, as the core- mantle interactions in our case, as part of the 
kinetic energy of Mercury. This section is widely inspired from iHenrardl (|200Sh . 

The components (vi, V2,v 3 ) of the velocity field at the location Xi inside the liquid core, in the frame of the principal axes 
of inertia of the mantle, are assumed to be: 



«i = \W2 + -V2JX3 - (w 3 + ^-uajxi, 

V2 = (u>3 + -2/3)2:1 — (<^1 + -^1^X3, 
V3 = (oJl + ^Uljx2 — (ui2 + —V2^Xl, 



(4) 

(5) 
(6) 



where (wi, W2, W3) are the components of the angular velocity of the mantle with respect to an inertial frame, and the vector 
of coordinates (ui, U2, v 3 ) specifies the velocity field of the core with respect to the moving mantle. 
The angular momentum of the core Nc is obtained by: 



TV^ = JJJ(^ x v)pdxi dx2 dx3 



(7) 



and the result is: 



We now set the following quantities: 

Di 
D 2 
D 3 

and we can write: 
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-v\ + wj I b + I -ui + wi ) c 



c \ 2 / « \ 2 

— V2 + UJ2 j a + I —U2 + UJ2 ] c 



b \ 2 f a \ ■ • 

-V3 + oj 3 la + I -V3 + lo 3 ] b 



fi 
h 



■^bc = 



(Ac -B c + C c ) (Ac + B c - Cc) 
-Ac + B c + Cc) (Ac + B C - C c ) 



■ab = M - Ac + Be + Cc) (Ac -Bc + Cc 



(8) 



iV^ = [A c ui + D lVl ] h + [BcU2 + D2V2] h 
while the angular momentum of the mantle is 



\CeUJ3 +D 3 U 3 ]f3 



(9) 
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N m = AmltJlfl + B m bJ 2 f2 + C m W3/3, 

and the total angular momentum of Mercury is 



(10) 



N = [Awi + D lVl ] /i + [Blj 2 + D 2 u 2 ] h + [Cuji + D 3 u 3 ] f 3 . 
The kinetic energy of the core is 



T c = i JJJ pv 2 dxi dx 2 dx s 



(11) 



(12) 



- (^A c (lu 2 + v\) + B c (lu 2 + f|) + C c (lu 2 + v%) + -Dicji^i + D2U2V2 + Dz^zv^j, 



while the kinetic energy of the mantle T m is 



-iV m ■ uj 



A 



From T = T m + T c we finally deduce the kinetic energy of Mercury: 



(13) 



(14) 



T = - (Aw 2 + Boj 2 + CujI + A c v\ + B c v\ + C c vi + 2LW 1 ^i + 2D 2 U2V 2 + 2D 3 oj a u :i ). 
We can easily check the expressions of the partial derivatives, as 

dT 



doji 



= Aji + Dii/i = Ni 



(15) 



(16) 



dT 



(17) 



where Ni are the components of the total angular momentum. N£ are not the components of the angular momentum of the 
core but are close to it for a cavity close to spherical. We have, for instance for the first component: 



N?-Ni" = (A B -Di)(ui-vi) 



— (c- b) (wi - vi), 
5 



(18) 



so the difference is of the second order in departure from the sphericity. From now on, we call angular momentum of the 
pseudo-core the vector iV c = ATf/J + N 2 fc + N£f 3 . 

With these notations, t he Poincare-Hough's equatio ns of motion , for t he system mantle-core in the absence of external 
torque, are (see e.g. Eq.15 in Touma and Wisdom (2001) or Henrard 1 2008! )'): 



dN _ 
IT = NXV * T ' 
dN c 



dt 
with 



Here T is the kinetic energy expressed in terms of the components of the vectors N and N c , i.e. 



(19) 
(20) 

(21) 



T=^- (A C N 2 + A{N{) 2 - 2DiNiN?) + ( B - N l + B ( N Sf - 2D 2 N 2 N%) 
Za Zp 

+-^(C C N% + C(N^) 2 - 2D 3 N 3 Ng) 



with a = AA C -D\,P = BB C - D\ and 7 = CC C - D%. 



(22) 
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Table 1. The shape parameters of Mercury. 



Parameter 


Value 


Reference 


J2 


(6.0 ± 2.0) X 10~ 5 


Anderson et al. (19871 


c 22 


(1.0 ±0.5) x 10~ 5 


Anderson et al. (1987) 


C/(MR 2 ) 


0.34 


Milani et al. (2001) 


5 = 1- C m /C 


0.421 


Mareot et al. (20071 


ei = J 2 MR 2 /C 


1.765 x 10~ 4 




e 2 = 2C 22 MR 2 /C 


5.882 x 10~ 5 





2.3 The Hamiltonian 

2.3.1 The rotational kinetic energy 

We assume that the cavity and Mercury are almost spherical, this allows us to introduce the four small parameters ti\ 

et = *L^Z*-J^, (23) 

e 2 = ^=2C7 22 ^, (24) 

2C C — A c — B c . . 

63 = 2C C ' (25) 

and also the parameter S — C c /C, i.e. the ratio between the polar inertial momentum of the core and of Mercury, ei represents 
the polar flattening of Mercury, while e 2 is its equatorial ellipticity. £3 and £4 have the same meaning for the cavity. If w e assume 



the co re of Mercury to be spherical, we should take £3 = £4 = 0, while £4=0 represents an axisymmetric cavity. iHenrard 



considered that the ellipsoid of inertia of the core and the mantle were aligned and proportional, the mathematical 



formulation was £3 = £1 and £4 = £2. Our parameters are g athered in Tab le [J 

We now introduce the two sets of Andoyer's variables I Andovei 19261 ). (I, g, h, L, G, H) and (l c ,g c , h c , L c , G c , H c ), related 



respectively to the whole Mercury and to its core. The angles (h,K,g) are the Euler angles of the vector n 2 , node of the 
equatorial plane over the plane perpendicular to the angular momentum N, the angles (J, I) position the axis of least inertia 
/1 with respect to 722- Correspondingly the angles (h c , K c , g c ) are the Euler angles of the vector n^, node of the equatorial 
plane over the plane perpendicular to the angular momentum of the pseudo-core N c , and (J c ,i c ) position the axis of least 
inertia with respect to n 2 . Figure [2] shows a schematic view of all the reference frames and relevant angles. The variables are 
(h, g, I) and (h c , g c , l c ) and the corresponding momenta (H — N cos K , G = N , L = N cos J) and (H c = N c cos K c , G c = N c , 
L c = N c cos J c ). Expressed in Andoyer's variables the components of N and JV C are: 

Ni = VG 2 - L 2 sin I, Nf = y/G 2 - L 2 sinZ c , 

N 2 = VG 2 - L 2 cos I, A" 2 C = yJG 2 -L 2 cos/ c , 

Nz = L, N§ = Lc. 

We can now straightforwardly derive the Hamiltonian Hi of the free rotation of Mercury, using Andoyer's variables and 
changing the sign of N c to take the minus sign of the Poincare-Hough equations into account (Eg I20fl . We also linearize the 
Hamiltonian with respect to the small parameters £; (their orders of magnitude being about 10~ 5 ), and get: 

Gl 



2C(1 - 


5)\ 


£1 




2C(1 - 


5) 2 


£2 




2C(1 - 


8) 2 


£3 




2C(1 - 


5) 2 


£4 





n o = T^r, K ( ° 2 + + 2 V(G 2 - L 2 )(G 2 - L 2 ) cos(l ~ l c ) + 2LL, 



^(G 2 - L 2 ) cos(20 + {G 2 C - L 2 .) cos(2; c ) + 2^J (G 2 - L 2 )(G 2 - L 2 ) cos(l + 1^ 



+2S^(G 2 - L 2 )(G 2 - L 2 ) cos(l + l c )J . (27) 



6 B. Noyelles, J. Dufey and A. Lemaitre 




Figure 2. The four reference frames gathered in the same view. The angles (h, K) position the plane orthogonal to the angular momentum 
N. The Euler angles (g, J, I) locate the axis of least inertia and the body frame (/i, /2, /j). The angles ( J c , l c ) place the angular momentum 
of the pseudo-core with respect to the axis of least inertia fi . 



We now introduce the following canonical change of variables, of multiplier n being the mean orbital motion of 
Mercury: 



p - 

r - 


= l + g + h, 
= -h, 


P = 
R = 


G 

nC ' 

P(l-cosif), 


$1 

Pc 

r c 


= - V /2P(1 - cos J) sinZ, 
= — l c + g c + h c , 
= -h c , 


Vi = 

Rc = 


v /2P(l-cosJ) cos/, 

G c 
nC ' 

= P c (l -COBKc), 


6 


= ^/2P C (1 + cos J c ) sin l c , 


V2 = 


^/2P C (1 + cos J c ) cos l c 



(28) 



In order to be consistent with the sign minus in the equations and before l c , the wobble of the pseudo-core J c has to be 
replaced by -k — J c . In this way, we have L c = G c cos(-7r — J c ) = —G c cos( J c ). In this new set of variables, we have 



m = nc^p?-(^~rf^, 
m = ncJpi-^^f^, 
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and the Hamiltonian of the free rotational motion becomes, after division by nC: 



+2(P 



Si /S2 +»?2 



)(&±*-*) 



nei 



2(1 - sy 



(&+nl pV , p2 ( p tl+ti 



+ V(p-^)(p,-^2)( M -e,6)) 



+ 



2(1 - sy\ I (4P " 5? ~ Vl) (c? " + 1 (4Pc " 6 " %) (6 " 772 } 



- 2V(p-^)(p.-^i)(^ + ^)) 



ne 3 



2(l-5) 2 



,2 £?+>J?\ 2 \ , / D 2 ,s1+% 2 



5 IP - (P- 



Pr ~ 



+ 25\[P 



+ 



(i ( 4P - £ - ^ ^ - + ( 2 - i ) I ( 4Pc - ^ 2 - ^ & - $ 



(29) 



Finally, in order to get an easy-to-use formula, we can develop this Hamiltonian up to the second order in (£1, £2, rji, 772) 



to get: 




■Hi = 


„ / 


2(1 - 5) { 


+ 


nei 


2(1-5)2 


+ 


nt2 


2(1-5)2 




ne 3 




2(1-5)2 


+ 


ne 4 


2(1-5)2 



{p 2 + ^f+ 2VPP c ( Vim - 6C2) + 2(7 



, 6 2 + ^2 2 | p jf + vl 



PPc 



p(£ + vi) + Pc{£ + m) + 2VpK{viV2 - 
J>(s1 - m) + p c (6 2 - - 2VPE( Vim + frfc) 

+ ^/i) + (2 - 7)^ (si + %) + 25VPP c (r, im - fcfc) 



+ (W, 2 - si) + (2 - i)p c (r, 2 2 - 6 2 ) + 2SVPPc(mm + £16)) . 



(30) 



2.3.2 The gravitational potential 

To compute the gravitational potential due to the Solar perturbation on Mercury we must first obtain the coordinates x, y, 
and z of the Sun in the reference frame linked to the principal axes of inertia (/1, /2, f 3 ). Five rotations are to be performed: 



/ Xi 



y j = R 3 (-l)R 1 (-J)R 3 {-g)R 1 (-K)R 3 {-h) 

z 



Vi 

Zi 



(31) 



with Xi, yi, Zi depending on the mean anomaly l , the longitude of the ascending node fl a , the longitude of the perihelion 
vj , the inclination i, and the eccentricity e. 
The rotation matrices are defined by 
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Rati) = 



COS ( 

sin <; 







/ 


1 








fll(0 = 







COS (j) 


— sin 




V 





sin <j> 


cos q 



The gravitational potential then reads: 

V 1 (l ,l,g,h,J,K) = 



(32) 



(33) 



where Q is the gravitational constant, M the mass of the Sun, (x,y,z) the unit vector pointing at the Sun in the frame 
(A) /si A) > while d is t he distance Sun -Mercury (expanded in eccentricity and mean anomaly). 

Let us note that unlike Henrard (2008), we consider that the perturbation is applied to the whole planet and not only to its 
mantle. We address the dynamical consequences later in the paper. 

From the variables x, y and z, it is easy to introduce the set of variables defined in (|28p . We also modify the moment A 
associated with l (that appears in the expressions of x and y) in such way that all our variables are now canonical with 
multiplier 1/nC and our gravitational potential becomes (after division by nC) 

H 2 {lo,P,r,R,ii,r]i) = ~\^^{ e ^ x2 +V 2 ) + ^(z 2 -y 2 ))- (34) 
Finally, we use the formulae (|30l) and (|34[) to get the Hamiltonian of the system: 



H = Hi (P, Ci, r,i , 6, m) + U-2 (lo, P, r, R,£i , Vi) • (35) 
The four degrees of freedom of this Hamiltonian are the spin (p, P), the obliquity (r, _R), the wobble of the whole body 
(£i, T]i) and the wobble of the core (£2, f)2)- In this study, we name "wobble" every motion dealing with a shift between the 
angular momentum of the body or its core, and its geometrical pole axis. It is different from the polar motion that concerns 
the rotation axis instead of the angular momentum. Contrary to the Chandler wobble for the Earth, we include in the term 
"wobble" every periodic contribution constituting this motion. 



3 COMPARISON BETWEEN AN ANALYTICAL AND A NUMERICAL STUDY 

To study this problem, we use both analytical and numerical methods that allow us to compare their efficiencies and check 
the reliability of the results. 



3.1 Analytical study 

In a previous paper by the authors IjDufev et al ] |2009l ). our model was a 2-degree of freedom Hamiltonian neglecting the wobble 
J, but including the planetary perturbations. Here we have a 4-degree of freedom Hamiltonian, but the way we perform our 
analytical study is similar to our previous paper. However there are some key diffe rences that we will highlight in this section. 
All the computations were made using our algebraic manipulator called MSNam 1 Henrardl 19861 ) . 



3.1.1 Resonant angles and Hamiltonian 



As mentioned earlier, it is a known fact that Mercury is in a 3:2 spin-orbit resonance. In other words, the rotation speed of 



Mercury p (where p = I + g + h, the spin angle of Mercury) is 1.5 times larger than its mean motion n , i.e. p ■ 



The angle 



la 



zu , with vj the longitude of the perihelion. The angle o\ actually represents the 



describing this resonance is o\ 
libration in longitude. 

The second resonant angle characteriz es the 1:1 commen s urabi lity between the orbital and rotational nodes, following the 
3rd of Cassini's laws ( Colombo ( 19661 ) or Lemaitre et al. ( 20061 ) for Cassini's laws applied to Mercury): 02 — r + fl a , with 
Q, being the longitude of the ascending node. This angle is actually linked to the latitudinal motion of Mercury (through its 
conjugated moment R which depends on the ecliptic obliquity K). 

Introducing the resonant angles in the Hamiltonian (|35p and using cartesian-like coordinates (expanded to order 5) for the 4 
degrees of freedom: 

x\ = o 1 (expanded around 0), y\ = P, 

X2 = y2P(l — cosK) sin <72, 2/2 
£1 = — y2P(l — cos J) sinZ, 771 
£2 = a/2P c (1 + cos J c ) sinZ c , 772 

the Hamiltonian is now H = H(l , x\, yi, X2, y%, £1, ?7i, £2, r]2). 

We also add constant precessions of the perihelion and the node, respectively dzu /dt = 0.2772831860533198 x 10 -4 rad/y and 



^/2P(1 — COS K) COS (T2 , 

V /2P(1 - cos J) cos/, 
y/2P c (1 + cos J c ) cos l c , 



(36) 
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d£l /dt = -0.2189047296429404 x 10~ 4 rad/y from the VSOP planetary theory ijFienga and Simonlliooj ) . These precessions, 



through the introduction of the resonant angle 02, will result in forced libration in latitude. 



3.1.2 Equilibria and free periods of the averaged quadratic Hamiltonian 

To compute the equilibria of the Hamiltonian, we first average the Hamiltonian over the fast angular variable (the mean 
anomaly Assuming that Mercury lies at the Cassini equilibrium and that there is no wobble motion (for the whole body 
and the core), we have xi — £2 = £1 = J?i = £2 = r\2 = 0. Putting that into the Hamiltonian, we compute the equilibria of 2/1 
and y2 using Hamilton's equations and an iterative process and we find y* = 1.5 — 6.117 x 10~ 7 and yj = 0.1502, resulting in 
an ecliptic obliquity of K* = 7° 1.873 arcmin. 

After a translation to this equilibrium, our quadratic averaged Hamiltonian looks like this: 

H 2 = ax\ + bx x x 2 + cx\ + dyx + eyiyi + fyl + g£i + h£i& + 1(2 + jr]i + kr]irj2 + lr)l- (37) 

We notice that the degrees of freedom related to the librations in longitude (xi,yi) and latitude (22,2/2), and those related 
to the wobbles of the planet (£1, r/i) and the core (£2, ^2) are coupled two by two, the coupling being much weaker in the first 
case than in the second one. 

To compute the fundamental periods of t his Hamiltonian, we must fir st disentangle the coupled degrees of freedom. To do 
this we use an untangling transformation ( Henrard and Lemaitrel 2005 ) twice, and after changes of variables to action-angle 
variables: 



xi = 

X2 = 
€l = 

e 2 = 



/2U sinu, 
/2V sinv, 
I2W smw, 
i2Z sin z, 



yi = v 2U cos u, 
2/2 = V2V cosv, 
771 = \J2W cosw, 
1)2 = s/2Z cos z, 



we have the following quadratic Hamiltonian: 



U2 = n u U + n v V + n w W + n z Z, 



(38) 



(39) 



with n u> n v , n w , n z the free frequencies corresponding respectively to the libration in longitude, the libration in latitude, the 
wobble of the planet and the wobble of the core. 

These frequencies (especially n v ) actually depend on the shape of the core. For the parameters given in Table [1] and £3 = ei 
and £4 = £2, the corresponding fundamental periods are 



T u = 12.0601 years, 
T w = 337.726 years, 



T v = 1065.99 years, 
T z = 58.6219 days. 



(40) 
(41) 



3.1.3 Lie perturbation theory 

To compute the evolution of the different variables, we use a perturbation theory by Lie transforms (see e.g. Deprit ( 19691 )). 



Our main Hamiltonian is the quadratic Hamiltonian described in the previous section and the perturbation contains all the 
other terms that we neglected to compute the fundamental frequencies (and to which we applied the same transformations). 
Here is a quick reminder on how this works. 
This perturbation theory is visualised through a Lie triangle: 

rto 

ill rto 



nil n/2 

rtl rto 

nil n_/3 

tt% Ttl rto 



(42) 



In the part Ho we put the averaged quadratic Hamiltonian explained in the previous section. In H\ we put what remains, i.e. 
the perturbation containing the higher order terms in U, V, W, Z, and the short-period terms. 

In the principal diagonal we will find the averaged Hamiltonian: H = X]"!'!) 01 ^oA- To get this averaged Hamiltonian we use 
the following homological equation: 

Ho = ht 1 + {K~ 1 \Wi) 

with the intermediate H" computed as follows: 



Tt . 



= H 



j+i 



i+E • )(n-l;Wi +i ), 
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where Wi is the generator of the ith floor of the Lie triangle, ( ; ) designates the Poisson bracket and I . ] is the binomial 



i 

coefficient. Note that the order (or the number of floors of the Lie triangle) is chosen in such a way that the transformation 
converges numerically, in other words we stop when we do not get more significant information by going one order further. 
These generators will help us to compute the evolution of any variable. Let us show how to get the generator of the first floor. 
We consider the first homological equation: Hq — Hi + (Hq;Wi). In this equation Hi is known and we choose Ho to be 
the average of Hi- In other words, Ho will contain only terms without short periods and of order larger or equal to 3. So, 
expanding the Poisson bracket, and using equation (|39p , we have the following equation to solve: 

dH° dWi &Hl dm | Q dH° dWi | Q dH° dWi | Q dH° dWi ^ ?/0 
du dU s^^y ^ u 'J^^Ls ^ v dw ^£Z^ dz o 1 

n u n v n w n z 

^ dWi dWi ewi dWi ! o 

<^ -ri u — n v — n w — n z — — = H ~Hi. (43) 

ou ov aw az 

Since Ho only consists of all the terms of H\ without short periods, the right-hand side term of the previous equation only 
contains short periodic terms. It is then easy to compute Wi and see that this is also only composed of short periodic terms. 
The computation of the other orders is done in a similar way. 

With the generators, we can now compute the evolution of any function of the variables. First, we go back to cartesian 
coordinates to avoid the singularities when any of the moment is (the angle is then undefined) and the formula is 

order 

f(xi,x2,£i,£,2,yi,y2,vi,m) = f(xi,x 3 ,ii,i2,yi,V2,m,ffii) + ~ C/C 215 !^ ^2, Ci, Ca, S/i, 2/2, T71, 772); Wi) , (44) 

1=1 

where the Poisson bracket is evaluated at the equilibria xi, £2, £i, £2, J7i, V2, fji, 772. 

After the use of this Lie algorithm, we have our transformed Hamiltonian in the diagonal, without short periods: H = 

E order -ui j -i 
i=0 U o/ %l 

Until here, except for the fact that we h ave 2 additional degrees of freedom, our process is very similar to the one described 
m our previous study (jPufev et al.ll2009h . The main difference is the change of fundamental frequencies. 



In this Hamiltonian, the linear terms in U, V, W and Z changed with the transformation process, yielding corrections to 
the free frequencies. These corrections also appeared in our 2-degree of freedom work, but so imperceptibly that we did not 
mention it. 

On the other hand, it plays a major role here. Here are the fundamental periods with the corrections: 

T u = 12.0568 years, T v = 1626.51 years, (45) 
T w = 337.853 years, T z = 58.6189 days. (46) 

We notice the very large change in T v , the period related to the libration in latitude. 

The period of rotation of Mercury is 58.646 days. The fundamental period related to the wobble of the core of Mercury is 
really close to this fundamental period and this particular combination of angles is present in the series related to the degree 
of freedom related to the libration in latitude. We are in fact close to a resonance that alters the efficiency of the algorithm. 
As a consequence, the numerical convergence of the algorithm is really slow when we take values of eg close to ex. For €3 = ex, 
we must use a 13th order Lie triangle to barely converge to the free period T v . During this process, we multiply series of 
millions of terms, which takes several days of computation. The process is divergent whenever £3 ^ 0.9 ex. 
Later in the paper, we draw a table of these periods for different values of £3 and £4 (cf. Table [2]). 



3.2 Numerical study 

As we already did for the rotation of Mercury with a spherical core jPufev et al.ll2009h or for rigid bodies (see e.g. iNovellesl 
1 20ld )L we performed a numerical study of the system by numerically integrating the equations derived from the Hamiltonian 
(|35[1 . and then carrying out a frequency analysis, for different values of the parameters £3 and £4. 

In order to integrate numerically the system, we first express the coordinates of the Sun (x,y) (as in Eal34[) thanks to 
Poisson series given by the VSOP planetary theory and the rotations given in (Ea l31[l , This way, we get coordinates depending 
on the canonical variables. Then we derive the equations coming from the Hamiltonian (|35[) . such as: 



dp 
dt 



dH 
dP' 



dP 
~dt 



dH 
dp ' 



(47) 



We then integrate over 13,000 years using the Adams-Bashforth-Moulton 10th order predictor-corrector integrator. 

The initial conditions are chosen close to the equilibrium and are iteratively refined to get amplitudes of the free librations 
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as small as possible (see e.g. Dufev et al. 1 20091 )). The reason is that we want to be able to extract the forced rotational motion 
of Mercury as accurately as possibl e, while the free terms are a source of noise. Moreover, these free terms are expected to 
be damped thanks to dissipations (|Pealel 2005). In addition to these forced terms, the proper frequencies are worth to be 
determined because of their significance on the dynamics of the system. We get their values from the first iteration. 

The frequency analysis algorithm that we use is based o n Laskar's original id ea, named NA FF as Numerical Analysis of 
the Fundamental Frequencies (see for instance Laskar 1 1993^ 1 for the method, and Laskar (2003) for the convergence proofs). 
It aims at identifying the coefficients a,k and uik of a complex signal f(t) obtained numerically over a finite time span [— T; T] 
and verifying 



/(*) ~ ak ex P(\/-lwfct), 



(48) 



where uik are real frequencies and complex coefficients. If the signal f(t) is real, its frequency spectrum is symmetric and 
the complex amplitudes associated with the frequencies LOk and —ujk are complex conjugates. The frequencies and amplitudes 
associated are found with an iterative scheme. To determine the first frequency U)i, one searches for the maximum of the 
amplitude of 



<j>{u) =< /(t),exp(v /Z Twt) >, 
where the scalar product < f(t),g(t) > is defined by 



(49) 



< f(t),g(t) >= ± j_J{t)g(t)* x {t)dt, 



(50) 



g(t)* being the complex conjugate of g(t) .x(t) is a weight function alike a Hann or a Hamming window, i.e. a positive function 
verifying 



2T 



x (t)dt = 1. 



(51) 



Using such a window can help the determination in reducing the amplitude of secondary minima in the transform (|50|) . Its 
use is optional. 

Once the first periodic term exp(v — lwit) is found, its complex amplitude a± is obtained by orthogonal projection, and the 
process is started again on the remainder fi(t) = f(t) — ai exp(V— luit). The algorithm stops when two detected frequencies 
are too close to each other, what alters their determinations, or when the number of detected terms reaches a limit set by 
the user. This algorithm is very efficient, except when two frequencies are too close to each other. In that case, the algorithm 
is not confident in its accuracy and stops. When the difference between two frequencies is larger than twice the frequency 
associated with the length of the total time interval, the d etermination of eac h fundamental frequency is not perturbed by 
the other ones. Although the iterative method suggested by IChampenois allows to reduce this distance, some troubles 

may remain. In our particular case, these problems are likely to arise because of the proximity between the free frequency of 
the core w z and the frequency of the spin. 



3.3 Results 



The Table [2] gives our analytical and numerical results, for 5 different sets of values for £3 and £4, ei, £2 being fixed with the 
known values of J2, C22 and 5=1 — C m jC. These 5 different sets are: 

• Case 1 : £3 = £4 = 0. This is a singular case, because the core-mantle interactions should not exist when the core is 
spherical. Moreover, we are at the exact resonance bet ween the proper mo de u) z and the spin frequency tj. We computed this 



case to look for an agreement with our previous study (|Dufev et al.| [2009), in which the spherical core was just removed. 
• Case 2 : £3/£i = 0.1, £4 = 0. We are close to the resonance, we here aim at detecting any discontinuity in the behaviour 

of the system close to the exact resonance. 

1. The cavity is homothetical to Mercury, this was the configuration studied by Henrard (2008) 



• Case 3 


£3/£l 


= £4/£2 


for Io. 






• Case 4 


: £3/£i 


= 3, £4 


• Case 5 


£3/£l 


= £4/£2 



3. By comparing this configuration with the previous one, we study the influence of the parameter 
£4, i.e. the equatorial ellipticity of the core. 



We made tests for a wider range of the parameters, but retained only these 5 cases for sake of conciseness. 
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Table 2. Proper periods of the system, numerically and analytically determined with a good agreement. Analytical values are missing 
when £3/ 1\ is small, because of singularities met by our algorithm of Lie transforms, uj is the spin frequency of Mercury, the last line 
represents the distance of the system with the resonance between u and ui z . We do not give any numerical value of T z for the exact 
resonance (£3 = £4 = 0) because we actually cannot numerically distinguish the free librations of the core from forced contributions at 
58.646 days. 
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Numerical 


Analytical 


Numerical 


Analytical 




Numerical 
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Numerical 


Tu (y) 


12.05800 


12.05775 


12.05685 


12.05772 


12.05685 




12.05777 


12.05685 




12.05773 


T v (y) 


615.77 


(large) 


1626.51 


1636.43 


1216.46 




1214.91 


1216.41 




1216.09 


T w (y) 


337.82 


337.82 


337.85 


337.87 


338.03 




338.14 


338.03 




338.20 


Tu (d) 




58.630 


58.619 


58.619 


58.585 




58.585 


58.585 




58.585 


Tz-uj (y) 




574.06 


344.88 


343.45 


154.08 




154.04 


154.05 




154.01 



Table 3. Evolution of the periods T v and T z with respect to £3/ei, with £4 = 0. These periods have been numerically determined, and 
confirmed analytically with a very good agreement when e^/ei > 1. 



C3/ei 


Tv (y) 


Tz (A) 


Tz-u (y) 


0.33 


3335.16 


58.628 


511.17 


0.7 


1966.31 


58.623 


409.08 


0.8 


1823.63 


58.622 


385.35 


0.9 


1718.34 


58.620 


363.50 


1.0 


1636.35 


58.619 


343.46 


1.1 


1570.86 


58.617 


325.10 


1.2 


1519.36 


58.616 


308.30 


1.5 


1408.10 


58.611 


266.01 


2 


1313.11 


58.602 


214.85 


2.5 


1250.26 


58.594 


179.64 


3 


1216.09 


58.585 


154.01 


3.5 


1198.68 


58.576 


134.72 


5 


1149.35 


58.550 


97.69 


10 


1107.62 


58.462 


50.83 



This table gives us first results on the influence of the shape of the core on the behaviour of the system. We can notice 
that the two periods T u and T m are quite constant. This is all the most interesting for T u because this degree of freedom, i.e. 
the longitudinal motion, is very weakly coupled with the variations of the obliquity, and even less with the two other ones. 
This means that this longitudinal motion is basically not influenced by the shape of the core, so the amplitude of longitudinal 
librations should depend only on C m /C, even with very accurate observations. 

The variations of the periods T v are worth noticing, because they present a discontinuous behaviour. From £3/61 = 3 to 
0.1 this period is getting larger and larger, some tests at £3/61 = 0.33 giving T v = 3, 335.16 years, while this period is too long 
to be numerically determined for £3/51 — 0.1. However, we have T v m 616 years at the exact resonance, what is consistent 
with the results obtained by simply removing the spherical core from the system. We think that it emphasizes the change of 
behaviour when the system is trapped into the resonance. As it can be seen in Table [3] and Figure [3] out of the resonance 
the period is decreasing, tending to reach the rigid value of 1,065 years. A least-square fit of T v gives T v « A{e^/e-i) B + C, 
with A = 564.488 ± 4.146 y, B = -1.25224 ± 6.003 x 10" 3 and ( 7= 1, 074.3 ± 3.233 y. By including ei in the constant and 
considering that 1,074 is close to the rigid value of 1,065 years ( D'Hoedt and Lemaitrel 2004 ; Rambaux and Boisll2004[ ). we 
can guess for T v an evolution such as 



T v ~ ^£3 ^ H~ T vri (52) 

where A is a constant and T vr is the rigid value of T v . The influence of £4 is very small. 

In the same way as T v , T z is decreasing when £3 is growing. It confirms that the resonance between the proper rotation of 
the core and the period of the spin is reached for £3 = 0. The third column of Table[3]is more significant because it represents 
the distance of this frequency ui z from the exact resonance. 
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Figure 3. Graphical representation of the proper period T v , associated with the obliquity. This curve can be fitted by f(es) ~ T vr + 
j4e 3 5 , where A is a constant and T vr is the value of T v when Mercury is considered as rigid. The circles are the outputs of the frequency 
analysis after numerical integration, while the solid line is /(es/ei) = A(e 3 / ei) s +C with A = 564.488 y, B = —1.25224 and £7 = 1, 074.3 

y- 

4 CONSEQUENCES ON THE OBSERVABLE ROTATION 

Our canonical variables are very convenient to describe the dynamics of the system, unfortunately they are not observable 
variables. Observations of the rotation of Mercury are in fact observations of its surface, i.e. the rigid mantle in our model. 
So, we have to express the components of the angular momentum of the mantle N m = A rn uifi + B m uj2f2 + C m uj 3 f3. 



4.1 The observable variables 

We deduce from the equations (|16[1 and (|17[1 : 



and 



DiNi + A°Ni 
D\ - AA C 



D 1 N 1 +AN? 
D\ - AAo 



(53) 



(54) 



and similar formulae for CJ2, 0J3, v 2 and vz- 

We can now easily deduce the expression of the angular momentum of the mantle with respect to the components of N 
and iV c : 



v . m DiW ~ A c Ni - pm D 2 m ^ B C N 2 - rim DsNS - C C N 3 r 

N = A ^5 — Jl+ti ^ =-5 J3- 



D\ - AA C 



Dl - BB C 



D\ - CCc 



(55) 



We can define a wobble J m and a precession angle l m of the mantle of Mercury this way: 



N m = G m sin J m sin; m /i + G m sin J m cos / m / 2 + G m cos J m f 3 , (56) 

where G m is the norm of N m . Because of the 3 : 2 spin-orbit resonance, G m is expected to be close to 3nC m /2. We get G m , 
J m and l m in equating the equations (|55[) and (|56|) . 

We now need to express the obliquity K m and the node h m of the mantle relatively to the inertial frame (ei,e2,e3). 
Naming Ti, T 2 and T3 the coordinates of N m in the inertial frame, we have: 
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t o 



= R 3 (h m )Ri(K m )R 3 (g n 



G" 



and we get: 

2i = G m s'm(Km) sm(h m ), 
T2 = —G m sm(K m )cos(hm), 
T 3 = G m cos(A' m ). 



(57) 



(58) 
(59) 
(60) 



Naming T[ c , T'2 and T 3 the coordinates of the angular momentum of the core N' c in the inertial frame (e"i,e2,e~3), we 



have: 



we now get: 



n c I = R 3 (h)Ri(K)R 3 (g)R 1 (J)R 3 (l) 



( K c 



(61) 



G m sin(A" m ) sin(/i„ 
G m sin(A ) cos (ft n 
G m cos(A 



and finally 

K m — arccos 
h m = arctan 



= Gsin(K)sm(h)-Ti c , 
= Gsin(K)cos(h) + T2, 
= Gcos(K) -Tz c , 



/ Gcos(A-)-T^ n 
V G m )' 
Gshi(K) sin(/i)-Tf 



(62) 
(63) 
(64) 



(65) 
(66) 



Gsin(A) cos(A) +T/, C 

From fEal57[). fEa !65p and fEa!66p. we can deduce g m . It is now straightforward to derive an angle p m , analogous to p, with 



Pm =lm+ 9m + h m - (67) 

The wobble J' c and the precession angle l' c of the core are directly derived from (Eq[9]), in the same way as J m and l m - They 
are not observable variables, but can be of planetological interest. 

4.2 Results 

The outputs given are chosen for their physical relevance. We express the longitudinal motion of Mercury, the obliquity of the 
mantle with respect to its orbital motion, the polar motion of the mantle and the wobble of the core. The first three variables 
can be observed, while the last one has indirect implications like on Mercury's magnetic field. 



4-2.1 Longitudinal motion 

The longitudinal librations have been studied for a spherical core jDufev et al.ll2009l : IPeale et al.ll2009h . These studies show 
two main contributions: a 88-day one of around 36 arcsec, due to the variations of the Sun-Mercury distance because of 
Mercury's eccentricity, and a 11.86-year one of approximately 42 arcsec, due to the Jovian perturbation on Mercury's orbit. 
Its amplitude is highly sensitive to the size of the core, because of the proximity of a secondary resonance with the proper 
frequency u u , the period associated being 12.06 years. We should keep in mind that the Jovian contribution has actually not 
been observed. Considering the uncertainty on the ratio (B — A) /Cm = (2.033 ± 0.114) x 10~ 4 , w e can only say that th e 
amplitude associated in the longitudinal motion of Mercury's mantle is bigger than 15 arcsec (see e.g. Yseboodt et al. ( 2010l ). 
Fig.5). 

We have determined the longitudinal librations thanks to a frequency analysis of the angle p m (Eq |67jl , after removal of 
a slope, i.e. the spin frequency of 39.1318408 rad/y. In all our numerical simulations, we get an amplitude between 35.82915 
and 35.830 arcsec for the 88-day contribution, and between 41.316 and 41.321 arcsec for the 11.86-year one. These results 
are in good agreement with the previous studies considering a spherical core, and do not show any significant variations. 
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Date (y) Date (y) 

E3/ei = ti/t'2 = 1 £3=3ei,<!4=0 

Figure 4. Variations of the orbital obliquity e of the mantle of Mercury for different shapes of the core. The time origin is J2000. 



Rambaux et al.l (|2007t ) had the same conclusion in a similar study, this means that the shape of the liquid (or molten) core 



cannot be derived from observations of the longitudinal motion. This result could be expected from the negligible variations 
of the proper period associated, i.e. T u (Table [2}. 



4-2.2 The obliquity of the mantle 



There are many ways to define the obliquity of the mantle. K m (Eq(SS| is the obliquity with respect of the inertial reference 
plane, i.e. the ecliptic at J2000. This quantity lacks of physical relevance, that is the reason why we here prefer to express the 
orbital obliquity e, i.e. the obliquity with respect to the normal to the orbit. It is derived from the scalar product between 
N m and the normal to the orbit, given by the cross product of the position and velocity vectors of Mercury. We show it in 
Figure H] 

This quantity shows essentially a secu l ar be haviour, as expected (see e.g. Peakj 1 20061) ; Yseboodt and Margot ( 20061 ) ; 
D'Hoedt and Lemaitrel d2008l ); |Pufev et all (|2009l )). Moreover, we can see from the plots that the obliquity is always close to 
1.66 arcmin, except in the strict resonant case (e-j = £4 = 0), where the obliquity is close to 1 arcmin. So, we can consider 
that there are two possible different behaviours: the resonant, and the non-resonant one. The resonant case is very improbable 
because it requires a strict fine tuning of the parameter £3 (related to the polar flattening of th e core). So, we can say that 
the shape of the core could probably not be derived from observations. As stated for instance in (|Pealdl2006l ) . these values of 
1 and 1.66 arcmin depend on the value of J2, on which we have a 30% uncertainty. 

We here consider the assumption of the liquid core to be valid at an y timescale. It is in fact often stated that Mercury 
should be considered as rigid, for the adiabatic evolution of the obliquity jPeale et al. 20021). Under this las t assumption, the 
expected obliquity should be about twice larger, as the 2.11 ± 0.1 arcmin measured by Mar got et al.l l|2007l ). 



4-2.3 Polar motion of the mantle 

We give here the polar motion of the rotation axis of the mantle about the geometrical North Pole. Following Henrardl (2005), 
we define the first two components Qi and Qi of the unit vector in the direction of the instantaneous axis of rotation by 
Qi Ri sin J m sinZ m [l + (C m — Am) /Cm] and Q2 ~ s in J m cos/ m [l + (C m — B m )/C m ], and multiply them by the polar radius 



of Mercury, i.e. 2,439.7 km (Sei delmann et al.| [2007) 
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Figure 5. Polar motion of the mantle, plotted over 5 years from J2000. 



Table 4. Synthetic representation of the polar motion of the mantle Q\ + y—XQ%. This motion can be expressed with a high accuracy 
(see the amplitudes) just with harmonics of the sidereal hermean frequency t. The two columns "Amplitude" arc related to the two cases 
shown in Figure [5] We can see that the differences between the two cases are quite small. 



T 


Amplitude (m) 
63/ei = £4,1^2 = 1 


Amplitude (m) 
£3 = 3ei,£4 = 


Period 
(d) 


1 


0.73344 


0.73118 


175.9 


-1 


0.44199 


0.44064 


-175.9 


3 


0.20077 


0.24330 


58.6 


-3 


0.17937 


0.17844 


-58.6 


5 


0.06283 


0.06264 


35.2 


-5 


0.06246 


0.06231 


-35.2 


7 


0.01765 


0.01760 


25.1 


-7 


0.01490 


0.01486 


-25.1 


9 


0.00507 


0.00505 


19.5 


-9 


0.00146 


0.00377 


-19.5 



In Figure [5] we show the polar motion of the mantle over 5 years, starting from J2000, for two different shapes of the 
core. In both cases, we see that this motion has a very small amplitude (smaller than 2 meters), and so should probably not be 
detected by the spacecrafts. They present similar aspects, the differences between the two representations being emphasized 
in Table H 

This table gives quasiperiodic representations of the polar motion of the mantle, defined by the quantity Qi + yf— TQ2, in 
the same cases as in Figure [5] The frequencies and the amplitudes associated have been numerically obtained using frequency 
analysis. The basic frequency of this motion, here named r, is the frequency of the sidereal hermean day. Its period is twice the 
orbital period of Mercury and thrice its rotational one. We can consider it as the frequency of precession of the rotation axis 
of the mantle about the geometrical pole. The other frequencies are harmonics of r. This table confirms that the amplitude 
of this motion is small. We can in particular notice that we cannot detect any resonant excitation of the 3r = cj contribution, 
while it is close to the free frequency u z . 



4-2.4 Polar motion of the core 

Even if the core cannot be observed, its rotation should still b e described. I t could indeed have some planetological implications, 
as for instance the origin of Mercury's magnetic field (see e.g. Christensen 1 20061 )). Figure [6] gives the evolution of its wobble J c 



for different values of the shape parameters £3 and £4, while Table [5] is a synthetic description of the quantity J c exp(y/—ll c ), 
representing the precessional motion of the core. 

The wobble J c is quite large when the system is close to the resonance, and its amplitude decreases when £3 (the polar 
flattening of the core) increases, i.e. when the system takes distance from the exact resonance. In fact, for £3 = 0, the behaviour 
of J c is a growing slope, meaning that the exact resonance tends to increase it dramatically, so the value of 40 arcmin that 
can be read on Figure[6]is not reliable. However, we are confident in the behaviour of the system out of the resonance, and we 
can see on the other three plots a quite similar visual aspect, except of the mean value. This aspect can be better understood 
thanks to Table [5] 

We see on this table the overwhelming predominance of the 58.6-d contribution, i.e. the rotational period. As expected, 
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Figure 6. Wobble of the core J c , obtained from numerical simulations after removal of the free librations. We can see that for a spherical 
core (t3 = £4 = 0) the visual aspect is very different from the other cases. Moreover, a long-term visualisation of J c shows a slope, i.e. a 
secular increase of the wobble of the core, while the other cases (out of the resonance) do not show it. 



Table 5. Synthetic representation of the precession of the rotation axis of the core about its geometrical pole axis. The quantity here 
analysed is J c exp(v / — ll c )- Contrary to the polar motion of the mantle, we see a high 58.6-d-contribution. The amplitude associated is 
raised by the proximity of the resonance with the proper frequency of the core. 



T 


Amplitude (arcsec) 


Amplitude (arcsec) 


Amplitude ( 


arcsec) 


Period 




e 3 = £i/10,£4 = 


£3Al = £4/c2 = 1 


£3 = 3<El 


e 4 = 


(d) 


3 


858.344 


85.532 




28.546 


58.6 


-3 




0.046 




0.010 


-58.6 


1 




0.049 




0.053 


175.9 


-1 


0.045 


0.011 




0.046 


-175.9 


5 


0.067 






0.005 


35.2 


-5 








0.002 


35.2 



it is excited by the proximity of the 1 : 1 secondary resonance with the proper frequency co z . We can see that the amplitude 
associated is roughly the mean value of J c as can be read from Figure [6] (858.344 arcsec = 14.306 arcmin). We also note some 
similarities with the precession of the rotation axis of the mantle (Table [4} , the frequencies involved being the same ones. 



5 CONCLUSION 

In this paper we have investigated the 4-degree of freedom behaviour of a rotating Mercury composed of a rigid mantle and 
a fluid ellipsoidal core, using both analytical and numerical tools with good agreement. We have emphasized the influence of 
the proximity of a resonance with the spin of Mercury, that can raise the velocity field of the fluid constituting the core. We 
cannot exclude a possibility of indirect detection of this effect by measuring Mercury's magnetic field. We have also shown the 
variations of the behaviour of the obliquity of the mantle with respect to the polar flattening of the core, this flattening being 
linked with the distance of the system from the resonance. These variations should be negligible unless the core is trapped 
into the resonance with the spin of Mercury. We have also shown that neither the observations of the longitudinal motion of 
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Mercury, nor of its polar one, could be inverted to get the shape of the core. However, they will give information on its size 
(i.e. the parameter S). 

Future works should take the viscosity of the fluid into account. It is assumed to alter the response of the core of Mercury 
to slow (i.e. « 10 J -year period) excitations, the planet being therefore assumed as rigid. As a consequence, a study of the 
spectral response of the rotation of Mercury on periodic solicitations with respect to the viscosity is worth studying. 
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